Phenomenological Theory of Superconductivity and Magnetism in Hoi_j;Dyj;Ni2B2C 
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The coexistence of the superconductivity and magnetism in the Hoi_i:Dya:Ni2B2C is studied by 
using Ginzburg-Landau theory. This alloy shows the coexistence and complex interplay of supercon- 
ducting and magnetic order. We propose a phenomenological model which includes two magnetic 
and two superconducting order parameters accounting for the multi-band structure of this material. 
We describe phenomenologically the magnetic fluctuations and order and demonstrate that they lead 
to anomalous behavior of the upper critical fleld. The doping dependence of T^ in Hoi_a;Dya;Ni2B2C 
showing a reentrance behavior are analyzed yielding a very good agreement with experimental data. 

PACS numbers: 74.20.De, 74.25.Dw, 74.25.Ha, 74.60., 74.62. 



The discovery of the series of nickelborocarbide, 
ENisBaC {R = Y, Yb, Dy, Ho, Er, Tm, and Tb) ||-|, 
a few years ago opened a new chapter in the discussion 
of magnetic superconductors |J]. This class of materi- 
als shows variety of magnetic ordered states |^ depend- 
ing on the rare-earth atoms R, because the magnetism 
is due to the spins of their localized 4/-orbital coupled 
via RKKY interaction. Some of these compounds show 
superconductivity with a comparatively high transition 
temperature {Tc — \Q.QK for LuNiaBaC). Band calcu- 
lations of non-magnetic LuNiaBaC pJlQl reveal a com- 
plicated multi-band structure with a rather high density 
of states at the Fermi level. Superconductivity is likely 
to be conventional originating from electron-phonon in- 
teraction due to the easily polarizable light ions B and 
C. The relatively large transition temperature is then a 
consequence of the enhanced density of states. 

An interesting aspects of this material class is that 
some members (i?=Tm, Er, Ho, and Dy) exhibit si- 
multaneously both magnetism and superconductivity. 
They have a rather wide range of the ratio of Tc/T^ 
from 0.6 for i?=Dy to 7.0 for R.=Tin. Among them, 
HoNiaBaC and DyNiaBaC constitute ideal materials to 
study the relation between superconductivity and mag- 
netism. They display the same antiferromagnetic or- 
der, but have different Tjv due to a different de Gennes 
factor {gj ~ lfJ{J + 1) (dGF). This leads to the sit- 
uation that the superconducting Tc lies above (below) 
Tjv for the Ho (Dy)-compound. Therefore, it is natu- 
ral to investigate how the two phases develop in the alloy 
Hoi_a:Dya;Ni2B2C flit]. The experiments show that there 
is no simple de Gennes scaling of T^ which is expected to 



monotonically decrease with increasing [gj — 1)^ J(J-f 1) 
according to the Abrikosov-Gor'kov theory, if we assume 
the R-ions enter as magnetic impurities. This is not the 
case here, since the R-ions form a regular lattice with 
magnetic order. In the crossing region of Tjv and T^ the 
superconducting phase boundary as function of doping 
concentration x is apparently discontinuous as a con- 
sequence of a reentrant behavior |ll| ]. The purpose of 
this letter is to explain these and further experimental 
features consistently within a single phenomenological 
model including both magnetism and superconductivity. 

The neutron scattering and magnetization experi- 
ments of HoNisBjC |l|] and DyNisBaC, show that the 
magnetic properties of the both materials are highly 
anisotropic with an easy axis in [110]-direction ^^■ 
Both HoNiaBaC and DyNiaBaC have antiferromagnetic 
(AF) order below Tm = 5 K and 10 K, respectively. 
The spins order ferromagnetically within the a-b plane 
and antiferromagnetically in c-direction. For HoNiaBaC, 
there exists a significant spiral short-range correlation 
(SP) with the wave vector, Q* ~ 0.91Q, where Q is the 
AF wave vector, i.e. the spins are aligned ferromagnet- 
ically in the a-h plane, and rotate by (/) ~ 163° between 
adjacent planes. A plausible microscopic origin of these 
two magnetic correlations has been discussed by Kalatsky 
et al. and Amici et al. Q. 

Based on this experimental information, we construct 
our phenomenological model of two magnetic order pa- 
rameters, Maf and Msp, corresponding to the two dom- 
inant correlations. The effective Ginzburg-Landau free 
energy has the form. 
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Tm = aAF{T)MAF^ + asp(r)Afsp' 
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where aAF{T) = aAFo(^-'?AF) and asp(T) = aspo(2"- 
Tsp). For the sake of simpUcity we restrict to the two 
dominant wave vectors (Q and Q*) of the spatial fluc- 
tuation only and take their mode-mode coupling into 
account |13]. We assume that the AF order is domi- 
nant (Tap > Tsp) and finally develops long-range or- 
der. Once the AF order is established, the SP order is 
suppressed, since the two magnetic orders compete with 
each other described by the mode-mode coupling term 
{Paf-sp > 0). Thus, only the AF order parameter has 
a non-zero mean value, as the temperature is lowered. 
Considering the fluctuation effects, we separate the or- 
der parameters into mean value and fluctuation part, 

Maf = Maf + SMaf and Msp = Msp + SMsp (2) 

where Msp — for all temperatures. The mean value of 
the AF order parameter is determined by minimizing the 
free energy with respect to Map, including the renormal- 
ization due to fluctuations. 



M 



aAp(T) + 3(3af{SMaf^) + Paf-sp{SMsp^) 



AP 



/?. 



AF 
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The Neel temperature Tjv is defined by the vanishing of 
Maf ■ In t^^ calculation of the mean square of the fluctu- 
ation in the transition region the renormalization due to 
the presence of the fourth-order (mode-mode coupling) 
terms has to be included. This can be done approxi- 
mately by applying a standard self-consistent decoupling 

6Maf^ = {SMaf^) + {SMaf^ - {6Maf^)), (4) 

SMsP^ = {SMsP^) + [SMsP^ - (SMsP^)), (5) 

leading to effective Gaussian fluctuation model [[l5[ . 
Thus, substituting Eq.(Q) and (^ into the free energy 
we obtain the following expression up to second order in 
the fluctuations of the order parameter. 



6Tm = («af(T) + 3/3afM1f + ^af{SMaf^) + I^afsp {SM sp^)) SMaf^ - ^{5M af') 



2\2 



(6) 



(asp(T) + pAF-SpMlp + /3sp(<5Msp') + Pafsp {SM af^)) 5Msp^ - ^(^Afsp')' 
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From this we can calculate self-consistently the Gaussian 
fluctuations including Eq.(|^) |15| 

(SMaf^) = fdSMAF SMAF^e"'^^^ / fdSMAFe-'^^^ (7) 

(<5Msp2) = fdSMsp SMsp'^e-P^^I fdSMsp e~^^^ (8) 

Since the modulation of the magnetization in the 
Hoi_:rDya;Ni2B2C occurs with a short length scale, i.e. 
has large wave vector Q or Q*, the superconductivity 
does not respond through coupling to the magnetic field 
(vector potential). For an (spin singlet) s-wave supercon- 
ductor the local spin polarization is pair breaking such 
that the above magnetic order should repel supercon- 
ductivity for electrons in orbitals subject to a net spin 
moment due to the 4/-spin ordering. (A quantitative es- 
timate of this effect requires rather detailed knowledge 
about the complicated band structure. Interestingly in 
experiment the Neel temperature does not deviate much 
from the de Gennes scaling, i.e. it seems to be basi- 
cally not affected by the superconducting order. Thus, 
we neglect the effects of the superconductivity on the 
magnetism.) Now we can calculate (A/^p) and (Afgp) 
using the above self-consistent scheme. The result of the 
temperature dependence of the intensities (M^) is shown 
in Fig.H. We choose the parameters to obtain a qualita- 
tive agreement with the neutron scattering experiments 
of the pure Ho-compound |q]. This yields a bare tran- 
sition temperatures Tap and Tgp which differ only by a 
few percent. Here the SP phase grows as temperature 
decreases and disappears quickly below the onset of AF 



order. 
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FIG. 1. Intensity of the antiferromagnetic A^af and c-axis 
spiral Afgp. The solid line denotes the intensity of the anti- 
ferromagnetic order and the dotted line the intensity of the 
spiral order. 

According to the band calculation of LuNi2B2C ||^, the 
Ni(3d) band has the largest contribution to the density of 
states at the Fermi surface suggesting that it dominates 
in the formation of the superconducting state. However, 
also the other bands involving hybridizations with B(2p), 
C(2p) and Ho(5d) or Dy(5d) should be included in the 
superconductivity. Hence, superconductivity of the nick- 
elborocarbide materials cannot be described by a single- 
band model. In connection with the magnetic order there 
is clear distinction among the bands depending on the lo- 
cation of the each element within the unit cell. For the 
AF order the influence of the rare-earth magnetic mo- 
ments is canceled on the Ni-sites located exactly in the 



center of a tetrahedron of the nearest Ho (Dy) atoms. 
The band originating from Ni(3d) does not feel the mag- 
netic momentum of Ho(Dy) anymore below the Neel tem- 
perature. On the other hand, the field generated by the 
moments of Ho(Dy) is not canceled at the site of the other 
elements. Moreover, because of large in-plane ferromag- 
netism of Ho(Dy) magnetic moments, the other ions feel 
a comparatively strong (pair breaking) spin polarization. 
In the case of SP correlation, the magnetic field at 
the Ni site is not canceled exactly. Thus, the SP phase 
would also affect the superconductivity of Ni(3d) band. 



For other ions the effect of the SP correlation is not so im- 
portant since the dominant effect comes eventually from 
AF order. Hence we propose simple phenomenological 
model which consists of two superconducting order pa- 
rameters. One superconducting order parameter (pA is 
associated with Ni(3d) band, and couples only to the SP 
order parameter. The other {^b) is due to the other 
bands involving essentially all elements. This supercon- 
ductivity is suppressed by the AF order. The Ginzburg- 
Landau free energy for the these two superconducting 
order parameter is 
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(9) 



where the parameters are chosen as Ti^a — 6.2K, and 
TcB = 2.9K (00 '■ standard flux quantum). 

For the superconductivity, we obtain two coupled dif- 
ferential equations for superconducting order parameters 



from the Ginzburg-Landau free energy (S) . If we assume 
Landau gauge, A — {0, Hx,0), then the linearization 
leads to the following homogeneous equations after solv- 
ing the differential equations. 
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The superconducting instability corresponds to the van- 
ishing of the determinant of these two equations which 
determines T^ and -ffc2- Note that in the absence of an 
external field the onset of superconductivity occurs at Tc 
higher than both Tca and Tcb as a result of the coupling 
of the order parameter components. 

Now we include the Dy-doping in HoNi2B2C. Dy and 
Ho has same magnetic ordering properties and simi- 
lar magnetic moments except for the magnitude of its 
dGF which determines the strength of the coupling be- 
tween localized 4/-spin and conduction electrons. Thus, 
we assume the doping of Dy changes only the average 
dGF which grows linear with the Dy concentration x 
(dGF(2;) = dGF(Ho) x (1 - x) + dGF(Dy) x x). We 
take 771 and 772 to be linearly proportional to the dGF. 
Note, that the Neel temperature due to RKKY interac- 
tion depends also linearly on the average dGF. 

As shown in Fig.g, our calculation explains well the 
experimental Tc variation even when Tc < T/v. At small 
doping {x < 0.2), the Tc determined by both the ipA and 
(fB, decreases due to the AF and SP fluctuations which 
also introduce an apparent discontinuity of the onset of 
superconductivity because of a reentrant normal state re- 
gion near the Tat. For x > 0.2 the upper superconduct- 
ing region ceases to exist. The lower Tc increases slightly 
with growing T/v and remains as the only superconduct- 



ing transition for larger Dy-doping concentrations. 
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FIG. 2. dGF and doping concentration vs. Tc for 
Hoi_iDyi:Ni2B2C. The dashed line denotes the Neel temper- 
ature and the solid hne the Tc. The dot-dashed line indicates 
the phase boundary in the reentrant region. The diamond 
points mark the Tc, the plus points the Tjv, and the filled 
diamond points the boundary of the reentrant region in the 
experiments [11]. 

Our two component order parameter model gains fur- 



ther support, if we compare the Hc2 curves of HoNi2B2C 
@ and DyNizBzC Q. FigJ| shows the Hc2 of our 
model for HoNi2B2C and DyNi2B2C. Our result repro- 
duces the dip in the upper critical field of Ho-compound 
near the Neel temperature as in Ref. flq]. Below the 
Neel temperature, both Hc2 curves are more or less iden- 
tical, because the remaining superconducting (Ni) band 
feels basically the same magnetic environment in both 
the Ho- and Dy-compound. 
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FIG. 3. Temperature vs. Hc2- The solid line denotes the 
Hc2 of HoNi2B2C, and the dotted line denotes the Hc2 of 
DyNiaBaC. 

As our model describes some of the physical proper- 
ties of the alloy Hoi_a;Dya;Ni2B2C very well, we may 
ask whether some conclusions can be drawn also for re- 
lated systems. Starting from DyNi2B2C, the Dy-ions 
may be replaced by other rare-earth elements. Since the 
superconductivity in this compound relies on the subtle 
canceling of the ordered magnetic moments on the Ni- 
site, other elements (different from Ho or Dy) disturb 
this balance substantially. Thus, doping DyNi2B2C with 
non-magnetic elements, Lu, (Lua;Dyi_:rNi2B2C) yields a 
large net moment on the Ni-site such that Lu should act 
like a magnetic impurity. Indeed the Tc decreases with 
the doping concentration of Lu in agreement with the 
Abrikosov-Gor'kov theory. In addition, Tm also drops 
due to dilution of the Dy-moments. A similar behavior 
is expected, if we dope the Dy-compound with a mag- 
netic element whose crystal field and spin-orbit coupling 
effects support a different ordering. Also such dopants 
disturb the balance and act as a magnetic impurity on 
the Ni-based superconductivity. Furthermore, these ele- 
ments introduce disorder in the 4f-spin system such that 
simultaneously T^v also decreases with doping concentra- 
tion. The above alloy gives results consistent with our 
assumption that the band relevant for superconductivity 
in the Dy-rich compound is Ni-based. 

From our model, which describes the behavior and 
mutual influence of two magnetic and two supercon- 
ducting order parameters, we can derive a consistent 
theory explaining all the key experiments in the alloy 
Hoi_2:Dy3;Ni2B2C. The multi-band structure is a crucial 
aspect to understand the superconducting properties in 
this system. It is an important feature that the Ni(3d) 



which dominates the superconductivity does not couple 
to the spin ordering of the rare earth ions. Ho and Dy. 
Thus, the mutual interaction of magnetic order and su- 
perconductivity is mainly due to the other bands which 
seem to be weaker superconductors readily suppressed 
by the magnetism. We have found that our model de- 
velops a reentrance behavior near the crossing region of 
Tc and T/y. This reentrance can also explain qualita- 
tively the non-monotonic temperature dependence of the 
electrical resistance for compounds in this doping region. 
It is basically a consequence of a sequence of normal- 
superconducting-normal-superconducting phase as tem- 
perature is lowered. Although the reentrant behavior 
was also reported in polycrystalline HoNi2B2C, our the- 
ory verifies that it appears without any inhomogeneity. 
Thus the reentrant behavior of Hoi^2;Dya;Ni2B2C is an 
intrinsic effect. This region may be very interesting to in- 
vestigate various fluctuation effects in and close to the su- 
perconducting phase. Clearly disorder effects must play 
an important role here too. Thus experimental stud- 
ies make only sense with high quality homogeneous sam- 
ples. In conclusion, we would like to emphasize that de- 
spite the fact that the superconducting phase is probably 
electron-phonon induced and conventional s-wave type, 
the physics resulting from the interplay with magnetism 
and the multi-band effects yield a wealth of unusual prop- 
erties reviving again the study of magnetic superconduc- 
tors. 
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